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ABSTRACT 

The possible existence of intermediate mass binary black holes (IMBBHs) in globular clusters (GCs) 
offers a unique geometry in which to detect space-time oscillations. For certain pulsar-IMBBH con- 
figurations possible within a GC, the usual far-field plane wave approximation for the IMBBH metric 
perturbation severely underestimates the magnitude of the induced pulsar pulse time-of-arrival (TOA) 
fluctuations. In this letter, the expected TOA fluctuations induced by an IMBBH lying close to the 
line-of-sight between a pulsar and the Earth are calculated for the first time. For an IMBBH consisting 
of IOMq and 10 3 M Q components, a 10 year orbital period, and located 0.1 lyr from the Earth-Pulsar 
line of sight, the induced pulsar timing residual amplitude will be of order 5 to 500 ns. 
Subject headings: pulsars :general — gravitational waves — black hole physics 
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1. INTRODUCTION 

The unique stability of the electromagnetic pulses 
emitted by radio pulsars could allow one to detect the 
presence of gravitational wave radiation. It is normally 
assumed that the pulsar, Earth, and all points along the 
line-of-sight between the pulsar and the Earth are lo- 
cated far enough away from the source that the space- 
time o scillations may be treated as gravitational plane 
waves (Sazhin 119781 iDetweilerl 119791. Recently, it has 
been suggested that intermediate mass black hole binary 
system s may form i n the centers of Globular Cluster s 
(GCs) IMUleJBOOa IWflll200l: iMiller fc Colbertl f20041. 
Given the high density of radio pulsars in GCs together 
with the relatively small scale length (100 pc) of a GC, 
it is possible to find a pulsar whose line-of-sight will pass 
closer to the binary system than its end points. It is 
even possible that the pulsar itself may lie close to the 
binary system. It is shown in this letter that this sce- 
nario differs from the plane gravitational wave case in 
two important ways. First, the amplitude of the induced 
residuals depends on the distance of closest approach, or 
impact parameter, between the binary system and the 
Earth-pulsar line-of-sight. Second, since the closest dis- 
tance may be well within one gravitational wavelength, 
the near-field terms will play the dominant role in induc- 
ing the observed pulsar pulse arrival time fluctuations. 

In this letter, the effect of the gravitational field of a 
binary system on pulsar timing residuals are calculated. 
This is the first time both near and far field terms are in- 
cluded in this calculation and no assumptions are made 
about the structure of the emitted gravitational wave 
front. The results are accurate to order (v/c) 2 and {d/r) 2 
where v is the characteristic relative velocity of the mass- 
energy within the source, c is the speed of light, d is the 
length scale of the binary system, and r is the smallest 
distance between the binary system and the Earth-pulsar 
line of sight. The general expressions for the timing resid- 
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uals derived here may be applied to pulsar timing data 
to constrain the properties of putative black hole binary 
systems. 

It should be noted that the 20 new millisecond pul- 
sars rece ntly discovered in th e dense globular cluster 
Terzan 5 (Ra nsom et al.ir2005|) may provide an interest- 
ing application of these results. The large number of 
pulsars together with this cluster's high central density 
(« lO 6 L0/pc 3 ) make it a good candidate to detect or at 
least limit the existence of an intermediate mass binary 
black hole (IMBBH) system in its core. 

In the next section, the general effect of a binary sys- 
tem on the arrival times of pulses emitted by a pulsar 
is calculated. We make use of geometric units where 
G = c = 1. Order of magnitude estimates of the in- 
duced timing residuals are also made for astrophysically 
relevant systems. The results are summarized in the last 
section. 

2. CALCULATING THE RESIDUALS 

The curvature of space-time will cause the observed 
pulsar frequency to vary as a fu nction of time. T he tim- 
ing residuals, R(t), are given bvlDctweilerl ll979() 



R(t) 



-dt 



(1) 



where t is time, v Q is the observed frequency, and v e is 
the emitted frequency. The frequency shifts will be cal- 
culated using perturbative methods. Given the four ve- 
locity of an observer, V 1 , the observed pulse frequency is 
given by v = —V^K^ where K v is the dual to the photon 
four velocity. Let V = W + 6V 1 and K M =K^ + SK^, 
where the barred quantities are the unperturbed values, 
and 6V and 8K are the corresponding perturbations. 
For the purposes of this discussion, we choose the un- 
perturbed four- velocities of both the emitter and the ob- 
server to be 1 in the time direction and zero for the three 
spatial components. To lowest order, the timing residu- 
als are given by 



R(t) 



6K$ 



dt . (2) 



/o K Q 
The last term on the right hand side is the standard 
Doppler shift due to the relative velocity between the 



2 



source and the emitter. The first two terms are deter- 
mined by the actual path of the photons as they travel 
towards the receiver. 

The geodesic equations will be used to solve for both 
5K and SV. To lowest order, these equations take the 
following form: 



dSK a 

dX 
d6V a 
dr 



1 



(3) 
(4) 



where A is the affinc parameter along the light-ray path, r 
is the proper time of the observer, and g^ v is the space- 
time metric. The metric is assumed to be of the form 
^jui/ + h^v where rj^ u is the flat Minkowsky metric and 
h^ v is a small perturbation given, at space-time location 
(*,**), by: 



(5) 



where T^ v is the trace-reversed stress energy tensor of 
the region generating the space-time perturbation. 

The timing residual due to ray path propagation will 
be calculated to first order in h using Eq. J2J rewritten 
in integral form as: 



sk§ - tig 



H(t + z,by + zz)dz 



(6) 



H(t, r) = -(/iqo.o + h„, + 2ho zfi ), (7) 

given in the black hole barycentered coordinate system 
shown in Fig. ^ with z pointing parallel to the light ray 
from the radio emitter to the observer (so that z = K°X), 
and x pointing from the black hole binary barycenter to 
the closest point on that ray. The emitter and observer 
positions z e and z a , and the impact parameter &, are 
defined in Fig. ^ 

Using Eq. JSJl, can be calculated using standard 
perturbative techniques to order (v/c) 2 and {d/r) 2 where 
v is the characteristic velocity of the source elements rel- 
ative to their center of mass, and d is the characteristic 
size scale of the source. Remarkably, H can be written 
to this order as H — dF/dz, where F is given by: 



F(t,r)- 



Qrr Q% Qrr Q zz ^Qy 



(r — z)r (r — z) 2 

+ Qzz - 2Q 

TZ 



(8) 



Qij is the second moment of the mass-energy distribution 
defined as 

(9) 



and the represents the time derivative. Given that 
the mass-energy distribution is usually specified in the 
(x, y, z) coordinate system, it is useful to note that 
Qrr = Qzz cos 2 6 + 2Q zx cos 6> sin 6> + Q xx sin 2 9, Q rz = 
Q zz cos 9 + Q zx sin6», and Q\ = Q xx + Q yy + Q zz , where 
r = V ' z 2 + b 2 and 9 = arccosz/r. 

Next, the contribution due to the 4- velocity of the ob- 
server and emitter will be calculated. Eq. Q determines 



the acceleration of the body in question. Under the as- 
sumption that the body does not move appreciably under 
the influence of the metric perturbation so that the time 
derivatives of r may be ignored, 



Kn 



G(t,r)-G(0,r) 



(10) 



where 



G(t,r) 



2Q r ,-(l + f)%-(l + f) 
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3Qrz ^>Qrr (l 



3Q r 



(i + ^)Q„ + (i + ^)f 



3Q rz 



Q\ 



15- 



Qrr 



dt' 



(11) 



Using Eqs. J2J), an( l (|10|l . the complete residual 
may be written as 

R(t) = /* {F(t', r Q ) - F(t let ,f e ) + 

G{t',r )-G(t Iet ,f e )}dt' (12) 

where t rc t = t' — [z — z e ] . Note that the terms which will 
ultimately give rise to secular terms (i.e. G(0, r)) are not 
included in the above expression. 

The above derivation assumes that the pulsar and ob- 
server are stationary in the barycentric frame of the bi- 
nary system. However, for relative velocities <C c, Eqs. 
|H1 ^2 an d ^] give the correct leading-order behavior of 
the residual, provided we allow z and r to be functions 
of t. 

For the case when z e goes to negative infinity and z a 
goes to positive infinity, one can show that the expected 
residual does not go to zero. Instead it limits to: 

Qxx(t ^o) Qyy(t To) 



R(t) = 2 



bit') 2 

Now when the source of the gravitational perturbation 
is a massive binary system, we can give an explicit form 
for the quadrapole moment in the barycentric frame: 

Qij = fis.Sj , (14) 

where s = s 2 — Si is the separation vector and = 
miiri2 / (mi + 7712) is the reduced mass of the binary sys- 
tem with component masses mi and 7772 . Using standard 
astromctric notation, we can write this directly in terms 
of the inclination i of the orbit to the plane of the sky, 
the position angle £1 of the ascending node, the argument 
of periapse u>, and the true anomaly v as shown in Fig. 
12 There r points towards the observer, <5 points to- 
wards the North celestial pole, s p is the direction of the 
separation vector at periapse, L points along the orbital 
angular momentum axis, and u points towards the as- 
cending node of the orbit. Defining a second unit vector 
w = L x u in the plane of the orbit, we can write the 
motion of the binary in that coordinate system as: 

s = a !+ecosv l cos (^ + v )u + sin (^ + V )M (15) 

where a is the orbital semi-major axis and e is the ec- 
centricity. Noting that x' — i cos 9 — z sin 9 is the pro- 
jection of x into the plane of the sky, and defining its 



dt' (13) 



position angle B in the same sense as fi, the separation 
vector in the (i, y, z) coordinate system is: 



the problem and that the residuals are proportional to 



I =*[ su cos(f! ■ 
+ y[ s u sin(f! - 
+ «[-,„ cos(S2 - 



5) cos e o - a m 
3) + s„ coa(n 
3) sin« + s w 
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Here we have taken v to be our independent parameter, 
which must be evaluated at appropriate retarded times. 
The time derivatives can be computed analytically using 
Kepler's law v = {2ir/P){\ + ecosu) 2 /(l — e 2 ) 3 / 2 , where 
P = 2-k^J a 3 /M t is the orbital period and M t = mi + m% 
the total mass of the system. Many standard techniques 
exist for computing the full time dependence of v(t), one 
of the more generic being to solve numerically the tran- 
scendental equation for the eccentric anomaly E: 



t-t p = ^-(E -esmE) 



(17) 



where t p is a time of periapse passage, and then to com- 
pute v via: 



tan(f) 



tan(f) 



(18) 



For simplicity, though, we will continue to write the ex- 
pressions in terms of (retarded) v. 

Solving explicitly for the case where b r e ,r a and b is 
effectively constant over the time of observation, Eq. I|13fl 
becomes: 



sin 2(u + v) sin 2(S7 — B) cos i + 



{sin 2 (w + v) cos 2 i — cos 2 (w + v)} cos 2(fi — B) 



(19) 



Next, order of magnitude estimates are made for the 
amplitude of the induced residuals using the above re- 
sults. Two cases will be considered. In case I, both the 
pulsar and the observer are infinitely far away from the 
gravitational wave source but the impact parameter is 
finite. In case II, the observer is infinitely far away, but 
the pulsar is at \z e \ ~ b <C P (i.e. the pulsar is in the 
near zone of the gravitational field). 

In case I, the amplitude of the induced residuals 
can be estimated using Eq. fTT^l : Ri ~ 2fia 2 /b 2 = 

2[aM^ 3 (P/2tt) 4 / 3 /b 2 . For the nominal case of a binary 
system with 10M Q and 1O 3 M components, a 10 year 
orbital period, and an impact parameter of 0.1 lyr, one 
obtains the following order of magnitude estimates: 



Ri ~ 5ns 



lOMe 



/ M t 
\ 1O 3 M 



2/3 



.1 lyr 



P \ 4/3 
^7r )<*> 



In case II, the residuals induced by the motion of the 
bina ry system will be dominated by the last term in 
Eq. i|ll|) . Hence an estimate for the residual amplitude is 

given by R H ~ (3/4)^M t 2/3 (P/27r) 10 / 3 /fe 4 . For the same 
system discussed above, one obtains the following: 



Rn ~ 500ns 



ioM Q y vio 3 M Q 



Mi 



2/3 



.1 lyr 



10 yr 



10/3 

(21) 



In order to understand the above scaling, note that the 
impact parameter, b, is the only external scale factor in 
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Fig. 1. — System configuration and coordinate system. 



the quadrapole moment Q ~ /ia 2 . A simple dimensional 
argument therefore gives R scaling as fia 2 /b 2 , times P/b 
for every time integral of Q in the leading-order term. 
Case I involves no time integrals while Case II has two. 
Kepler's law is then used to write a in terms of the orbital 
period. 

3. APPLICATION AND DISCUSSION 

General expressions for the periodic timing residuals 
induced by a binary system were calculated. It was 
shown that systematic variations in the pulsar timing 
residuals depend not only on the location of the pulsar 
and the observer, but also on how close the binary sys- 
tem is to the pulsar observer line-of-sight. As long as 
the line-of-sight impact parameter is finite, a non-zero 
residual amplitude can still occur even if both the pulsar 
and the observer are infinitely far away from the binary 
system. For a given impact parameter, the residuals cal- 
culated using case I and case II represent the range of 
possible residual amplitudes provided that z e < 0. In 
other words, the pulsar must be near or behind the bi- 
nary system. 

Globular clusters present an interesting opportunity 
to discover intermediate mass binary black holes using 
pulsar timing. Table ^ shows the timing residuals, for 
cases I and II, that a 10M Q + 10 3 M Q binary system with 
a 10 year period would induce on known pulsars near 
(as specified by the impact parameter b) to the cores of 
their respective clusters. All pulsars in the table have 
periods < 10 ms. By way of comparison, the root-mean- 
square (RMS) timing noise for millise cond pulsars is ap- 
proac hing the level of 100ns or better l)van Straten et all 
|2001|) . New efforts like the Parkes Pulsar Timing Array 
(PPTA) project are actively working to improve the RMS 
noise level 4 . The proposed square kilometer array (SKA) 
project will provide timing precisions as low as 10 ns in 
the next 10-20 years and will also uncov er all pulsars 
beamed at Earth within globular clusters l)Cordes et alJ 
2004; IKramer et aLll200l . 

Of course, a single pulsar could never definitely detect 
a binary system, although it could be suggestive. In or- 
der to make a strong case, timing residual oscillations 
must be seen in two or more pulsars and these oscilla- 
tion must be consistent with the same binary system. 
The 20 milli second pulsars in Ter zan 5 may offer such an 
opportunitv ljRansom 

Part of this research was carried out at the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der a contract with the National Aeronautics and Space 
Administration and funded through the internal Re- 
search and Technology Development program. 
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TABLE 1 

Residuals induced in globular cluster pulsars by a IOMq + 10 3 Mq, 10- year binary 



Globular Cluster 


Pulsar 


6(lyr) 


Rj (ns) 




47 Tuc 


J0024-7204O 


0.26 


0.8 


11 


47 Tuc 


J0024-7204W 


0.34 


0.5 


4 


NGC 6266 


J1701-3006B 


0.18 


1.6 


50 


NGC 6624 


B1820-30A 


0.37 


0.1 


3 


M28 (NGC 6626) 


B1821-24 


0.12 


4 


200 


NGC 6752 


J1910-5959B 


0.38 


0.4 


2 


NGC 6752 


J1910-5959E 


0.49 


0.2 


0.9 


M15 (NGC 7078) 


B2127+11D 


0.19 


1.1 


10 


M15 (NGC 7078) 


B2127+11H 


0.37 


0.4 


3 



R/I (ns) 



Reference 



Freire et al. f 20031 
Camilo et al. (20001 
Possenti et al. (20031 
B|'jl.'> H al i ! 
Rutledge^^l^ (2004) 
D'Amico et al. (2002) 
D'Amico et*aL f2002i 
Anderson ( 19931 
Anderson Q9931 




Fig. 2. — Binary orbital parameters. 
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